GAUSSIAN HYPERGEOMETRIC SERIES AND SUPERCONGRUENCES 



ROBERT OSBURN AND CARSTEN SCHNEIDER 

Abstract. Let p be an odd prime. In 1984, Greene introduced the notion of hypergeometric 
functions over finite fields. Special values of these functions have been of interest as they are 
related to the number of F p points on algebraic varieties and to Fourier coefficients of modular 
forms. In this paper, we explicitly determine these functions modulo higher powers of p and 
discuss an application to supercongruences. This application uses two non-trivial generalized 
Harmonic sum identities discovered using the computer summation package Sigma. We illustrate 
the usage of Sigma in the discovery and proof of these two identities. 



1. Introduction 

In [16] and [T7], Greene defined general hypergeometric series over finite fields. His aim was 
to show that these functions satisfy properties analogous to classical hypergeometric series. For 
example, the four major evaluations of the ordinary hypergeometric series 3-F2 due to Saalschiitz, 
Dixon, Watson, and Whipple [5] all have finite field interpretations (see page 126 of |16j). 
Greene's work was in part motivated by the analogy between Gauss sums and Gamma functions 

BD1, M, M- 

His approach has proven to be a powerful technique for character sum evaluations. Recently, 
several authors have shown that special values of these functions are related to the number 
of points over ¥ p , p an odd prime, of Calabi-Yau threefolds [3], traces of Hecke operators |llj . 
formulas for Ramanujan's r-function [35], and the number of points on a family of elliptic curves 
|13j . We should also mention that hypergeometric series over arbitrary fields has been developed 
[H], [E], but their application to number theory has yet to be investigated. 

The purpose of this paper is to further study arithmetic properties of hypergeometric functions 
over finite fields. In particular, we explicitly determine these functions modulo higher powers of 
p and then briefly discuss extensions of supercongruences. 

We first recall some definitions. Let ¥ p denote the finite field with p elements. We extend all 
characters \ of F* to ¥ p by setting x(0) := 0. Following [16] and [17], we give two definitions. The 
first definition is the finite field analogue of the binomial coefficient. If A and B are characters 
of F p , then 
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where J(x> VO denotes the Jacobi sum if x and ip are characters of ¥ p . The second definition is 
the finite field analogue of ordinary hyper geometric functions. If Aq, A\, . . . , A n , and B\, . . . , 
B n are characters of ¥ p , then the Gaussian hyper geometric function over ¥ p is defined by 

« t :::: 

where the summation is over all characters x of F p . In this paper, we restrict our attention to the 
case Ai = <p p for all i and Bj = e p for all j where <j) p is the quadratic character and e p is the trivial 
character mod p. We shall denote this value by n+ \F n {\). By [IB] and [T7], p n n+ iF n {\) 6 Z. 

Before stating the main result, we recall that for i, n E N, generalized Harmonic sums -£/n are 



defined by 

n 



7 

i=i J 



and fl™ := 0. We now define the quantities 



X(p,A,n) :=0 P (A)^(^ +J ) (Y) (-1)^ ( 1 + 2(n + " ff f 

(3) j=o V J / V / V 2 

+ J - 2 (^(l + n)(^ + . -Ff ) 2 - (^)(H^ + .-H?))\ 

— i i 

Y{p,\n):=^ p {\)j2( E ^ +J ) (*f) (-iy l \-™(l + (n + l)j(HW 

(4) j=° V V 2 +J 



2 /n \ 21 



(5) Z{p, A, n) := ^ p {\) £ ( 2j ) 16^' A"^ 



and 



P-5 

j' 2 



(6) *>(■>■*):=£ , i+M 2 » + 1 » 2A ' J "- 

j=0 IIi=0r 2^ 



where Z := ^4-^. The main result of this paper is the following. 
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Theorem 1.1. Ifn>2, then 

(7) -p n n+1 F n {\) = (-M-W +1 [p 2 X(p,\,n) +pY(p,X,n) + Z(p,A,n)] (mod p 3 ) 
and if n = 1, t/ien 

(8) -p 2 i ? i(A) =p 2 [x(p,X,l) + D(p,X)j + pY(p,\,l) + Z(p,\,l) (modp 3 ). 

We note that Theorem 11.11 generalizes both Theorem 1 in p] , where the case n = 2 was 
handled modulo p 2 , and Theorem 2.4 in [25J. As an application of Theorem II. 1[ we prove a 
supercongruence for the Legendre symbol (-^)- This result generalizes Theorem 1 in [27j . 

Corollary 1.2. Let p 6e an odd prime. Then 

The method of proof for Theorem 11.11 has its origin in [3] . Namely, the idea is to first observe 
that since the functions n +iF n (\) are defined in terms of Jacobi sums, then one can express 
them as Gauss sums. One then applies the Gross-Koblitz formula [18] to express the Gauss 
sums in terms of p-adic Gamma functions. Using combinatorial properties of the p-adic Gamma 
function, Theorem 1 1 . 1 1 then follows. For an introduction to these methods, see |32| . This general 
framework has been the basis for several recent results on supercongruences (see [I], |21j . [25], 
[27J, [28J, |29|). Theorem 1 1 . 1 1 has recently been used to settle a conjecture of van Hamme (see 
[26]). Finally, the congruence in ([9]) appears to hold modulo p 4 . This has been numerically 
confirmed for all primes less than 5000. 

The paper is organized as follows. In Section 2, we recall properties of the p-adic Gamma 
function. In Section 3, we prove Theorem ll.il In Section 4, we prove Corollary II .2l using Theorem 
11.11 and two non-trivial Harmonic sum identities discovered using the computer summation 
program Sigma [39]. A description of the non-trivial methods involved using the Sigma package 
is included in Section 5. We should also mention that similar harmonic number identities were 
discovered and proven in [33J . These types of identities played an important role in the proof of 
Beukers' supercongruence for Apery numbers (see [2] or Theorem 7 in [3]). 

2. Preliminaries 

We first recall the definition of the p-adic Gamma function and list some of its main properties. 
For more details, see [22], [30], or [37]. Let | • | denote the p-adic absolute value on Q p . For 
n e N, we define 



r p (n):=(-l)» II 

j<n 

CmO=i 

One can extend this function to all x € Z p upon setting 



-1 
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T„(x) := lim T„(n). 

n—>x 

The following Proposition provides some of the main properties of T p . 
Proposition 2.1. Let n G N and x G Z p . Then 

(i) r p (o) = i. 

/o% r p( x + !) f -» i/|x| = 1, 



v ; T p (x) " 1 -1 if\x\ < 1. 

(3) 7/0 < n <p - 1, then n\ = (-l) n+1 r p (n + 1). 

(4) |r p (x)| = i. 

(5) Let xq E [1, 2, . . . ,p] be the constant term in the p-adic expansion of x. Then 

r p (x)r p (i-x) = (-ir°. 

(6) If x = y (mod p n ), then T p (x) = T p (y) (mod p n ). 
For x G Z p , we define 



(10) Gi(x) 
and 

(11) G 2 (x) 



r p (x) 
r p (z) 



r p (x)' 



One can check that C?i(x) and (^2(2) are defined for all x G Z p using the fact that T p (x) is 
locally analytic and |r p (x)| = 1. We now mention some congruence properties of the p-adic 
Gamma function. For a proof of this result, see [8] or |21j . 

Proposition 2.2. Let p > 7 be prime, x E Z p , and 2 G pZ p . TTien 

(1) G a (x), G 2 (x) G Z p . 

(2) We We 

T p (x + z) = T p (x) ^1 + zGi(x) + yG 2 (x)^ (mod p 3 ). 

(3) T' p (x + z) = F' p {x) + zF' p f (x) (mod p 2 ). 

We also need the following combinatorial congruence which relates T p to certain binomial 
coefficients. 

Lemma 2.3. If p is an odd prime and 1 < j < , then 
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Proof. By Proposition 12.11 (3) and (5), we have 



J! 2 (^-J)! 

r p (i+j + 2)r p (i+j-|) 



r P (i+j) s 



Now, using Proposition 12.21 (2). we have 



{r P (I + ,) + |r;(I + J )}{r P (I + ,)-|r;(I +J )} (mod 



p 2 ) 



r p(^+/) (modp 2 ) 



and the result follows. 



Finally, we need to define 



□ 



(12) A(j):=Gi(£+i)-Gi(l+i) 
and for a positive integer n 

B(n,j) :=^(G 2 (i + j) - G 2 (l + j)) + + jf 

+ ( " +1) 2 (w+2) Gi(l +3? ~ (n + l) 2 Gx(i + j)Gi(l + j). 

We require the following Lemma which relates A(j) and B(n,j) to generalized Harmonic 
sums. The proof is similar to Lemma 4.1 in [21J and thus is omitted. 

Lemma 2.4. Xe£ p be an odd prime and < j < 2=-. T/ien 



r=0 

and 

(15) B(n,i) - +i - Hf) 2 - m(H^ + . - Hf) (mod P ). 
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3. Proof of Theorem 11.11 
We are now in a position to prove Theorem ll.il 
Proof. Let n > 3 be odd. From ([I]) and ([2]), we know that 



^+ A(A) = -^^J(0,x^ +1 x(A) 

1 — V 



p 

x 

where \ is the complex conjugate of x- After expressing the Jacobi sum J(<j), x) n+1 in terms of 
Gauss sums, we then apply the Gross-Koblitz formula [18] to get (see also [3] or [25j l 



co ] (X) 



- P \ + iF n (X) = — I P (A) + (-0 P (-1))— £ - -f^j 

(16) 1 V -° PU 

+ n+i y ^ (A) I 

. p+1 L P\p-l 2^ /J 
J ~ 2 

s — 1 

Here u; is the Teichmiiller character which satisfies u>(X) = X p (mod p s ) and thus 

J (A) = X jpS ^ (mod/) 

for s > 1. As n > 3 is odd, the second sum in (|16j) vanishes modulo p 3 . As j^— = 1 + p + p 2 

(mod p 3 ) and thus = — j — jp — jp 2 (mod p 3 ), we apply parts (5) and (6) of Proposition [2J] 
and reindex the summation to obtain 

(17) 



p\ +1 F n (X) 



p-1 



- (i + p + p 2 ) P (A) + (-^ P (-i))— £ rnl-I I ( mod A 

[ r p (i + j + jp + 2P l ) n+v j 

By Proposition 12.21 (2), we see that 

r p (x + j + jp + jp 2 ) n+1 = T p (x + j) n+1 [l + (n + l)(jp + ip 2 )Ga(z Q + j) 

+ ^O'P + JP 2 ) 2 (g 2 (x + j) + ndixo + j) 2 )l (mod p 3 ) 



for xo € Z p . We expand the numerator and denominator of (|17p with xo = \ and xo = 1 
respectively. After multiplying the numerator and denominator by 



1 - (n + l)jpG 1 (l+j) - 2±1jY(g 2 (1 +j) - (n + 2)Gi(l + j) 2 ) - {n + l)jp 2 G x {\ + j), 
we get 
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p n n+1 F n (X) 



(18) 



2-1 

i2 ±i ^r p (i+j)" +1 



(1 + P + fc(A) + (-fc(-l)) V £ _^Z^_(i + ( n + l) iM (j) 



p-l 



+ (n + l)jp 2 ^(j) + j 2 p 2 B(n,j)juj 2 J(A)| (mod p 3 ) 
where A(j) and B(n,j) are denned by (|12p and (|13p . We now need to consider the sums 



(19) P (A) + (-^(-l))— £ / ; + J n+1 (l + 2(n + l)jA(j) +i 2 B(n,i))a; 2 ^(A) ; 



^ T p (l + j) 



2-1 

n+ l 2 p (I _|_ 7 N ) n + 1 / 

(20) 4> P w + (-^(-i))— e r p ; ( ; + jj n+1 (i + (n + 

and 



l)iA(j))o;V-^A), 



p-i 



(21) 



n+l 2 p /l + -)n+l ! 

>(A) + (-0p(-i))— E /nJU ^' 3 ^ 



which are the coefficients of p 2 , p, and 1 respectively in (|18p . Observe that as we want to 
determine n+ iF n {\) mod p 3 , it suffices to compute (fT9|) mod p, ([20]) mod p 2 , and ([2Tj) mod p 3 . 
Also note that 



(22) 



wV-J'(A) = wV(A)ur-»'(A) = (j) p (\)uj-i (X) 



as a; is of order p — 1. By Lemma 4.4 in [21], we see that 



j'-i 



(23) " s " 2 ( n + WpE Ta^TF (mod p2 

By Lemma |2.3| we have 



r=0 



(24) 



U-l)(-l) 



n+l 



r P (|+ir +1 
r p (i+j) 



n+l 



(mod p 2 



and so after combining Lemma 12.41 (|22p . (|23p. (|24p and accounting for j = 0, then (|19p is 
congruent modulo p to ([3]) and (|20l) is congruent modulo p 2 to ([4]). Here we have used the fact 
that T p (l) 2 = 1 and T p (\) 2 = -</> p (-l) and thus for n > 3 odd 
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TJ±) n+1 n+l 



r P (i) 

By Proposition 2.5 in [27], we have 

and thus using (|2"2"j) . we have that (|2"Tj) and ([5]) are congruent modulo p 3 , namely 

P-i 

p (a) + (-0 P (-i))— E T^irf^"" (A) 
( 25 ) i=1 " ^ 

2 /o \ n+l 

3\ 



3=0 

This proves the result for n > 3 odd. A similar argument applies to the case n > 2 is even. 
We now turn to the case n = 1. By (|16p . we need only consider the last sum 



(26) -M-W E rfff/ M- 

. p+1 i pU-i 2^ 



By (5) and (6) of Proposition 12.11 and after reindexing the exponent of uj(X), (1261) is equivalent 
modulo p 3 to 

By repeated use of Proposition 12.11 (2) we have for 2+1 < j < p - 2 that 



1 p+l\2 

2 + ~ 

1 p + 1 \2 _ /p + 2\2 



r P U + ^H =1 (modp) 



Mh^^) 2 -(^) 2 (^) 2 



r,(i + ,- 2 )' = (^)'(^)'...(E±i)" 

By Proposition O (3), T p (l + j) 2 = (j\) 2 . Also using the fact that X p 1 = 1 (mod p) , we have 
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3p— 3 3p— 3 

uj 2 -J(A) = A 2 -i (mod 



p-i 

EE A~ 



-p-i-j 



(mod p) 



P-l , 
= A 2 (mod p) 



for 2±1 < j < p — 2 and thus 



P ^ ^3_ 



J - 2 



' ^ + i)V 2 (A) + ... + ^ 



(2+1)12 V2 



p+2\2 /2p-5\2 



(p-2)P 



2 p+1 



_2_w -+p-2j w 2 (A) (modp) 



1 



(1) 2 (2)2... (E^)2 



,(-!)(*) 



^(-l)(p-2)2(p-3)2...(2±l)2V 2 



p-3\ 2 1 P-5 

A 2 (mod p) 



= -<^(-l)D(p,X) (modp). 
This proves the result for n = 1. □ 

4. Proof of Corollary 11.21 

Theorem 11.11 can be used to obtain modulo p 3 super congruences in various settings. For 
example, Apery numbers [3], traces of Frobenius endomorphisms on elliptic curves [24J, [31] . 
and colored partition functions [31] all occur as special values of n +i-F n (A) for certain n and A. 
We do not mention these results here, choosing instead to illustrate with one example. We now 
prove Corollary 11.21 

Proof. If p is an odd prime, then by Section 3 in [16] , 





X(p,l,l)+D(p,l) +pY(p,l,l)+Z(p,l,l) (mod/). 



We now claim that 



(27) X(p, 1, 1) + D(p, 1) + 1 ee (modp). 

In order to verify ([27]) . we first study X(p, 1, 1). By ©, we have 
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X{P , m) = g ( % + J ) (T) c-tf ( x + - 



The identity 



(28) g(-l) fe + fe ) Q (l + 2k{H n+k ~H k ))= (-l)»(2n + 1) 



was discovered using Sigma (see Lemma 2.2 in |27j). In particular, we find 
p-i 



i=o 

and thus 

(29) 

X(p,l,l) 



p-i 



+ E ( % + (f >-^ 2 (K*^ - ^') 2 - - "?>))) (-d ,). 

By Proposition 12.21 and Lemma 12.31 we also have 

P-3 _ x _ x 

(30) D(p,l) = E^ + ^ (V) ( mod ^)- 

For positive integers n, the relation 



(31) 



+e("; j )"(")" ( - i), =" ( - i+2 " )( - i) " 



was found using Sigma. Equation ([27]) now follows upon taking n = 2^- in (f3Tj) . in the identity 
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02) e( --i+c-iy 



n + j\ ( n N 

and reducing modulo p. We now consider Y(p, 1,1). By we have 




+j 3 



•./(//; . // ; 



For positive integers n, the relation 



(33) 



t("; , ')G)(- i ) i ( i ^-Ki + <- 2ff i i) )) 

- v l + 2n) J (-1)" - -n(l + n)(-l)" £ ^ 

^ ' i=l 

was discovered using Sigma. Taking n = in (|3"3"j) and reducing mod p 2 , we have 

£-1 

(34) Ffe i ) i)= p+ 3 ( _ 1) ^^W1 (modp2) _ 

i=l ^ ' 

Equation © then follows from ©, (JSTJ) , and (|MD- □ 

5. Finding and proving identities (f5Tj) and ([33)1 with Sigma 

An efficient algorithm to find and prove identities involving nested definite and indefinite sum 
expressions, such as (|3ip and (|33p . has only recently been developed and implemented. Inspired 
by hyper geometric summation [53], in particular Zeilberger's creative telescoping method |42j . 
and Karr's indefinite summation algorithm |19|, 120] (extended to definite summations), the sec- 
ond author developed and implemented an algorithm using Mathematica to handle various 
summations. The resulting package is called Sigma. For a more detailed description of the 
algorithms incorporated into Sigma, please see [39]. Applications of this computer algebra pack- 
age include proving identities that arise in the enumeration of rhombus tilings of a symmetric 
hexagon [121 138] . in the verification of Stembridge's totally symmetric plane partitions theo- 
rem |41] [6], and in certain Pade approximations [9j. In this section, we illustrate how the 
package can be used to discover and prove identities (f5Tj) and ([55]) . For simplicity, we write 
for 
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5.1. Identity (I33|) . With Sigma we produce the following harmonic sum identities: 



(35) ( n + Q kH k =(-l)"n(n + l)(2H n - 1), 

(36) f>l) fc ( H + k ) (fj kH n+k ={-l) n n{n + l)2H n - (-l) n n 2 , 
and 

(37) f>l)* (" + k k ) (t) kH - k = ~ + ^ + C" 1 )^ 2 " + ^ 
fc=o ^ / V / 



2n N 

n /2i\ 



o n /2i\ 

+ 2n(n + l)(-l) n H n - -n(n + !)(-!)» ^ ^ 



i=l 

Then, combining ([32]) . (|36|) and ([37|) we arrive at identity ([33]) . 



Remark 5.1. Note that (|28p . (|35p and (|36p can be proved using hyper geometric techniques 
which appear in [5j and [7j. The key observation is that differentiation of the rising factorial 
(x)k = x(x + 1) . . . (x + k — 1) (resp. l/(cc)fc) in x and afterwards substituting x = 1 produces 
(l)k-fZfc (resp. —Hk/(l)k). With this fact, one can produce, e.g., ([35]) by setting up the identity 

(-n) fc (n+ l) fc n(n + 1) n(n + 1) (a? - n - l) n -i 
( 38 ) L fci737 fc = — 2 i ; i(l-n,n + 2;x + l;l) = 

with Gauss' theorem, differentiating (|38|) in x, and setting x = 1. Obviously, the successful 
application of this technique relies on the fact that one knows the underlying hypergeometric 
identity such as (|38|) for the particular case (|35p . It would be interesting to see proofs of 
identity (|37p . in particular of identity (|3ip . along the lines sketched above. Recently, a skillful 
application of partial fraction decomposition has been used in [36] to derive identity (|37|) . but 
so far no proof of identity (|3ip has been found. 

Subsequently, we illustrate the computation steps for identity (|371) which can be executed in 
a straightforward manner. After loading the package 

ln[l]:= << Sigma. m 

Sigma - A summation package by Carsten Schneider © RISC-Linz 

into the computer algebra system Mathematica, we insert the sum in question: 

ln[2]:= S = SigmaSum[SigmaPower[— 1, k]kSigmaBinomial[n + k, k] 

SigmaBinomial[n, k]SigmaHNumber[n — k], {k, 0, n}] 



k 



Remark 5.2. Various functions support the user, like SigmaSum for sums, SigmaPower for 
powers, SigmaBinomial for binomials, or SigmaHMumber for harmonic numbers. 
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Next, we compute a recurrence relation for the given sum S by inputting: 
ln[3]:= rec = GenerateRecurrence[S] 

0ut[3]= (n + 2)(2n+ l)(n+ l) 2 SUM[n + 2] + 2(n + 3)(2n 2 + 4n + l)(n + l)SUM[n + 1] + 

(n+ l)(n + 2)(n + 3)(2n + 3)SUM[n] == -(2n+ l)(2n + 3)(3n+ l)(3n + 4)(-l) n 2n 

n 

This means that SUM[n](= S = J2(-l) k k( n + k )(l)H n _ k ) satisfies 0ut[3]. 

k=0 

Proof of 0utf3j: Define f(n, k) := (— l) k k( n ~t ) The correctness follows by the creative 

telescoping equation 

(39) g(n, k + 1) - g(n, k) = c (n)/(n, k) + ci(n)/(n + 1, k) + c 2 {n)f{n + 2, fe) 

and the proof cer tificate c (n) = (n + 2)(n + 3)(2n + 3), ci(n) = 2(n + 3) (2n 2 + 4n + l), 
c 2 (n) = (n + l)(n + 2)(2n + 1) and 

flf(n, fc) = (fc - l)k 2 (2H n _ k (k - n - 2)(fc - n - l)(ra + l)(fc(4rc + 7) - 2(2n 3 + 10n 2 + 17n+ 10))+ 
(-k -n- l)(16n 4 + 88n 3 + 179n 2 + Win + 2k 2 (4n 2 + Yin + 7) - k(24n 3 + 

98n 2 + 131n + 59) + 58)) (-l) fe + ^ Q / ((n + l)(-fc + n + l) 2 (-fc + n + 2) 2 ) 

delivered by Sigma. We verify f|39|) as follows. Express g(n,k + 1) in terms of h(n,k) = 
(-l) k ( n t h ) (£) and # n _ fc by using the relations 

(n - k)(n + k + 1) 
(fc+1) 

and 

H n -k-l = H n -k r • 

n — re 

Similarly, express /(n + i, k) in terms of /i(n, k) and H n -k by using the relations 

+ 1, k) = H — 7 — /i(n, ft) 
n — re + 1 

and 



^ fc + 1) = - v ,7; ; ^ (n, *o 



n — A; + 1 

Then (i39l) can be checked directly. Summing (i39|) over k from to n produces Out[3]. □ 

Next, we solve the recurrence relation Out[3] by typing in: 
ln[4]:= recSol = SolveRecurrence[rec[[l]], SUM[n]] 



n 1 

Out[4]= {{0, n(l + n)(-l) n }, {0, (1 + n)(-l) n ( - 1 + 2n^ -)}, 

i=i 

{1, 4(-l) n ( - 2(1 + 2n) M + 3n(l + n) ]T ^)}} 
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The result has to be interpreted as follows. Sigma finds two linearly independent solutions 

n 1 

h\(n) = n(l + n)(— l) n and ti2(n) = (1 + n)(— l) n ( — 1 + 2n^^— ) of the the homogeneous 

i=l 1 

version of 0ut[3] (indicated by the in front) plus one particular solution 



1 /2n\ 

p(n) = - 2 (! + 2n ) ( n J + 3n ( 1 + n ) E 



i=l * 

of the input recurrence itself (indicated by the 1 in front). The correctness of the result can be 
easily verified by using, e.g., the relation 

( 2i \ 1 _ ( 2i \ 1 2 ( 2n + X ) ( 2n 
^W^~^W^ + (n+l) 2 \n 

Finally, by taking all linear combinations c\hi{n) + c<ih<i{n) + p(n) for constants c\ and C2, 
free of n, we obtain all solutions of Out[3]- Hence, by considering the first two initial values of S 
we can discover and prove ([3T]) : 
ln[5]:= FindLinearCombination[recSol, S, 2] 

Out[5]= -(-l) n (l + n) 2 + (-i)'(l + 2n) M + 2n(l + n)(-l) n £ j - |n(l + n)(-i)» £ ^ 



Remark 5.3. Looking at the identities (|32p . (|36p and (|35p one immediately sees that the com- 
bination ([32]) + 2 • (f36|) — 2 • ([35]) produces ([28]) . Since the sums can be combined so nicely, we 
had also the hope to find a solution for the sum 

S x (n) := ±(-X) k + & ) Q (l + 2k(H n+k - H k 

Sigma was able to compute the recurrence 

(n + 2) 2 S x (n) + (2A - 1) (in 2 + 18n + 21) S x (n + 1) 

+ (l6n 2 A 2 + 80nA 2 + 100A 2 - 16n 2 A - 80nA - 100A + 6n 2 + 30n + 39) S x (n + 2) 

+ (2A - 1) (in 2 + 22n + 31) S x (n + 3) + (n + 3) 2 S x (n + 4) = 0, 

but failed to find any solution for a generic value A. Interesting enough, choosing A = \ the 
recurrences gets much simpler. In particular, this indicates that considering the sums S\(2n) 

and Si(2n + 1) separately, one can compute recurrences of order 2 for each of them. Indeed, 

applying the mechanism from above for each of the sums gives (two different) recurrences of 
order two. Luckily, we can even solve the recurrences which yields 

2 V ; (2n)\ 

and 

S \ (2n + 1} = ( 2^(ff ! ( (2n + 1} {Hn ~ H2n) ~ l ) ■ 
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5.2. Identity (j3Tj) . Finally, we derive the two identities 

£(" t *) (*) (_1>i ( 2k2 ( ff " +t " " fc2 (" s ' " fff ' , 

(40) =(1 + n) 2 (-2 - 2n + n 2 ) 1 / 1 {, + -n(-4 + lira + 6n 2 + 3n 3 )(-l) n 

2(2 + 2n)\ 4 

" ^(-1 + n + n 2 ) + \n\l + „)'(-!)" £ + n 2 (l + n) 2 (-l)" £ 

i=i ^ i=l 

and 

-i / \ -i 



(41) - ((1 + n) 2 (-2 - 2n + n 2 )\^_J_ + _ n (_ 4 + Un + 6n 2 + 3n 3 )(-l) n 

1 3 n -|2 n (_-\ 

" + n + n 2 ) + ^n 2 (l + n) 2 (-i r £ T^yT + ^ + ») 2 (-l) n £ ^ 

i=i ^ ' i=i 



which immediately gives identity (|3ip . 

One option is to follow the same strategy as above: We can compute a recurrence for 

,n [6 ] : = mySum = |>l) k ( n + k ) M ( 2 k 2 ( Hn+k - H k ) 2 - k 2 ( H £ k - H< 2 > 

and can solve the derived recurrence to find the right hand side of (00]). But, since the found 
recurrence relation is rather big (it has order 4), and the proof certificate is even bigger (it fills 
about one page), we follow a refined strategy presented in [33] and [40] . Namely, by running 
our creative telescoping algorithm with the additional option Simplif yByExt — ► DepthNumber 
we can find a recurrence of smaller order (order one!): 

In[7]:= rec = GenerateRecurrence[mySum, SimplifyByExt — > DepthNumber] 

Out[7]= 2(2n + l)(n + 2) 2 SUM[n] + 2(2n + l)n 2 SUM[n + 1] == 

4(1 + 2n) + n 2 (n + l)(n + 2)(3n + 2) ]T ^llM® + 2n (4n 2 + 3n - 4) (2n + 1) ^(-l) 1 ® + 



8(n - i)n(n + l)(n + 2)(2n + 1) ^-l) 1 ^ 1 ) 0H n+i - ^(-i) 1 ^ 1 ) ©Hi 



Proof of OutUS: Define f(n, k) = (-l) fc ( n f ) (£) (2k 2 (H n+k - H k f - k 2 {H^ k - flf )). Then 
the correctness of Out[7] follows by the creative telescoping equation 



(42) 



g(n, k + 1) - g(n, k) = c (n)f(n, k) + ci{n)f{n + 1, k) 
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with the proof certificate co(n) = 2(n + 2) 2 (2n + 1), c\(n) = 2n 2 (2n + 1) and 

g(n, k) = U(k - l) 2 n(n + l) 2 (2n + l)k 2 (Wf - AH n+k H k + 2H 2 n+k + ff< 2 > - H™ k 

(n(n + 2)k z - (n 3 + 2n +2n + 2)k 2 - (n + l) 2 {n 2 - 2)k + n(n + l) 2 (n 2 +n- 2)^8n(n + 1) 
(2n + 1) ( H k ~ H n+k ) + (16n 5 + 48n 4 + 29n 3 + 14n 2 + 20n + 8)fe 2 + n(n + l) 2 (16n 4 + 23n 3 + 



n + Yin + 8) - (32n b + 101n s + 98n 4 + 55n a + 54n 2 + 36n + 8)k ) , \ fc Kk} N + 

1 (1 - k + n)n(l + n) 



n 2 (n + l)(n + 2)(3n + 2) £ \^ )a AJ + 2n (4n 2 + 3n - 4) (2n + 1) (?) 

i=0 j=0 

k k 

+ 8(n - l)n(n + l)(n + 2)(2n + 1) (^(^TCV) (?)#»+* ~ E^TO C) H *) 

i—0 i=0 

Since the sums and products inside of g(n, k) are all indefinite, e.g., we can apply the relation 

fc+l k 

Et-tfPrt (?) = Ec-tfrfX?) - (-if^fflO. 

i=0 i=0 

the verification of ()42p is immediate. Summing (|42p over A; from to n produces Out[7]. □ 

At first glance the recurrence 0ut[7] seems to be disappointing: we start with the definite sum 
mySum, and end up with a recurrence again involving definite sums. But, these sums are much 
simpler than the input sum. In particular, facilitating again Sigma, we can produce mechanically 
the identities 



n (n+i\ fn\ 



nl 2 



(n + i) 2 n 2 {2n) 



i=0 
and 

g-)' (" : *) (i) * - d-d* c : *) (i) - 

Using in addition (|32p . we can simplify the recurrence to 

ln[8]:= rec = rec/.{f:(-l) i mP - -(-1)--^, ^(-l) 1 ^) (?) - (-1)", 

i=0 i=0 

E(-i)'m(?)H. - E(-i)'m(?)Hn+.} 



(-l) n (n+l)(n+2)(3n + 2)(n!) 2 „ , 

4(2n + 1) - — A ^ + 2(-l) n(2n + 1) (4n + 3n - 4) 



Out[8]= 2(2n + l)n 2 SUM[n + 1] + 2(n + 2) 2 (2n + l)SUM[n] == 

[n + 2 
(2n)! 

Given this recurrence, one can directly read off its solution. With some simplifications Sigma 
yields: 

ln[9]:= recSol = SolveRecurrence[rec[[l]], SUM[n], SimpleSumRepresentation — > True] 
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Out[9]= {{0, (-l)V(n + I) 2 , {1, (1 + n) 2 (-2 - 2n + n 2 )||^| - ±(-1 + n + n 2 ) + 



in(-4 + lln + 6n 2 + 3n 3 )(-ir + \ n 2 (l + n) 2 (-l)» £ T^y + + ^f^T £ ^}} 

i=l ^ i=l 

Looking at the first initial values we end up at the identity (|40p . 

Finally, we attack the sum S(n) on the left hand side of (|4ip . Namely, we derive the recurrence 

( ±o\?( ^ (-l)"(n + l) 2 (n + 2)(3n + 2)(n!) 2 

(ra + 2)5(n) +n S'(n) = - — 2; 

we remark that for this hyper geometric sum any implementation of Zeilberger's algorithm 
could do the job. Using this information, we find as above, the right hand side of (|4ip . 
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